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I. INTRODUCTION 



The mystery of "flavor" poses a major challenge in particle physics. The unpredicted masses 
and mixings of three families of leptons and quarks in the flavor sector compose thirteen out of 
nineteen free parameters in the whole Standard Model (SM). Weak scale supersymmetry (SUSY) 
[|l|] as a leading candidate for new physics beyond the SM provides no further understanding about 
the origin of flavor. In fact, it extends the mystery of flavor by necessarily adding three families of 
squarks and sleptons. Without additional assumptions for flavor structure of the soft SUSY break- 
ing, supersymmetric theories often encounter phenomenolgical difficulties, known as the SUSY 
flavor problem 0. In dynamical models of electroweak symmetry breaking, the phenomenolog- 
ical constraints on flavor-changing-neutral-currents (FCNC) make it hard for model building to 
accommodate the observed heavy top quark mass, unless new dynamics associated with the top 
quark is introduced || . Predicting the full mass spectrum and mixing pattern in the flavor sector 
may have to invoke new physics scales ranging from the weak scale up to very high scales in a 
single unified theory. Less ambitious and more practical approaches follow a "bottom-up" path, 
which effectively parametrize the new physics with flavor in a way not explicitly invoking unknown 
dynamics at the high scales. For instance, certain ansatze for fermion masses matrices were ad- 
vocated ||], and realizations of horizontal symmetries were proposed^, || to explain the fermion 
mass hierarchies and Cabibbo-Kobayashi-Maskawa (CKM) mixings @ at relatively low scales. The 
FCNC fermionic Yukawa couplings to the Higgs sector can also be constructed in a phenomeno- 
logically viable way || |9| . When the electroweak symmetry breaking sector is nonlinearly realized 
flTo| , |TT| , U ], new FCNC gauge interactions of fermions can be economically described by effective 



operators of dimension-4 [13, 14]; while for the linearly realized Higgs sector, these couplings arise 



from dimension-6 effective operators [15|. 



The recent exciting evidence for neutrino oscillations 1 16 , 17] strongly points to nonzero neutrino 
masses with large (rather than small) mixings which further discriminate the lepton sector from 
the quark sector and deepens the mystery of "flavor". In fact, the atmospheric oscillation data Ji~6[ ] 
favors maximal mixing between the /j, and r neutrinos [^l] via the Maki-Nakagawa-Sakata (MNS) 
matrix [18] [52]. This development has led to a great amount of theoretical effort, with the hope of 
revealing the underlying new physics in the leptonic flavor sector [|Hi|]. 

In this paper, we systematically explore the low energy constraints on effective operators induced 
by large r— fi mixing. This is strongly motivated by the neutrino oscillation data, in particular, the 
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favored maximal mixing between the second and third generations via MNS matrix. Note that the 
large or maximal mixings in the MNS matrix may come from either neutrino mass diagonalization, 
or lepton mass diagonalization, or from both of them. It is indeed tempting to search for charged 
lepton flavor violations, in addition to the existing neutrino oscillation experiments. There are two 
classes of structures that can lead to flavor mixing for charged leptons. The first class is that there 
is no tree-level flavor mixing for charged leptons after mass diagonalization. This is an analogue 
to the SM quark sector, where the flavor mixing effects appear at tree level only in the charged 
current sector via the CKM mixing matrix. For instance, in a SM-like framework with additional 
mass and mixing parameters in the neutrino sector, r — fi mixing is generated at one- loop level and 
is generally suppressed by a factor of m^/M^, which would be negligible. The second class yields 
tree-level mixing effects after the mass diagonalization from the flavor eigenbasis into the mass 
eigenbasis. Typical theories in this class include extended models with extra Z' or Higgs doublets, 
generic weak-scale SUSY models, and dynamical models with compositeness, which often have rich 
structures of flavor mixing, leading to testable new phenomena. By contemplating on the large 
lepton mass hierarchy, dimensional analysis suggests the new physics scale associated with the 
mass and mixing of the third family leptons to be the lowest one in the lepton sector. 

Before experiments can directly access the new physics scale associated with the lepton flavor 
sector, we use an effective theory formulation, obtained by integrating out the heavy degrees 
of freedom from a more fundamental theory. In Sec. II, we present our general formalism for 
the dimension-6 fermionic effective operators involving r — \i flavor mixing with typical Lorentz 
structures 

Gurr)(rrY), (i.i) 

where T contains possible Dirac 7-matrices. This is beyond the simplest flavor-diagonal form of 
the effective four- Fermi contact interactions pQ| . With the operators in Eq. ([L"l|), we analyze their 
relationship to the realizations of the SM gauge symmetry and the underlying flavor dynamics. We 
further estimate the expected size of their coefficients, for a given cutoff scale at which the effective 
theory breaks down and new physics sets in. We also comment on to what extent our formalism 
can be applicable to loop-induced processes. In Sec. Ill, we systematically explore the constraints 
on the new physics scale associated with each operator, mostly from current experimental bounds 
on rare decays of r, hadrons or heavy quarks. For the operators involving at least one light quark 
(u,d,s), the current low energy constraints typically push the new physics scale to a few TeV or 
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higher. Those operators with two heavy quarks are at the present subject to weaker constraints, 
only about a few hundred GeV. Some of the scalar and pseudo-scalar operators are free from 
any experimental constraint so far. We summarize our results in Table I and conclude in Sec. IV. 
Appendix A analyzes the relationship between our operators introduced in Sec. II and those with 
explicit SM gauge symmetry in both linear and non-linear realizations. In Appendix B, we further 
extend the above Eq. ( |0| ) to include the forms with generic lepton bilinear £ a T I 13 , and derive the 
corresponding constraints from rare r decays. 

II. FORMALISM FOR EFFECTIVE OPERATORS WITH r-/x FLAVOR VIOLATION 
A. Constructing the Dimension-6 T — fx Operators 

We consider an effective theory below the new physics scale A, which can be generally defined 

as 

Aff = £sm + AC (II.2) 

where Csm is the SM Lagrangian density and AC denotes the new physics contribution via effective 
operators. For the current study, we will focus on the dimension-6 operators involving third and 
second family leptons (r, fi), 

AC = A45 = J2^-(t r, r) (r T 3 gP) + H.c. , (II.3) 

where Tj € (1, 7 5 , j a , 7^75) denotes relevant Dirac matrices, specifying scalar, pseudoscalar, vector 
and axial vector couplings, respectively. We will not consider the possibility that T has tensor 
structure due to the following reasons: (a) there are no two-body decays involving the tau with 
tensor structure, and thus any bound would be extremely weak; (b) the tensor matrix elements 
involving quark bilinears are either unknown or known very poorly, and thus bounds would be not 
only very weak, but also rather uncertain; (c) the tensor structure does not generically appear in 
most models that we know. Here, we consider Tj to be the same for both t—/j, and q a — q@ bilinears, 
where a and (3 run over all allowed combinations of quark flavors. Under these considerations, we 
can show that Eq. ( |II.3| ) is the most general form (containing one r — fj, bilinear and one quark- 
bilinear) which respects the unbroken C/(l) em gauge symmetry. As shown in the AppendixA, 
the operator Eq. flH.3|) corresponds to a nonlinear realization of the electroweak gauge symmetry 
under which all fields feel only the unbroken U(l) em . For the linearly realized electroweak gauge 
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symmetry where the physical Higgs and would-be Goldstone bosons form the usual Higgs doublet, 
the dimension-6 r — fi operators have to respect the full electroweak gauge group Gsm = SU (2) l <8> 
U{\)y and are thus found to have a more restricted form, 

A£ li„ea r = A£ Uneax = £ _g_ ( _ p . ^ ^ p . ^ + jj c {RA) 

where the chirality indices x> C = L, R. This restricts T £ (-y a , 7 CT 7 5 ), so that Eq. ( |II.4| ) only belongs 
to a sub-set of operators in Eq. (II. 3). [53] At first sight, this is somewhat surprising as the scalar and 
pseudoscalar operators in Eq. ( |II.3|) are fully absent at dimension-6. However, it is interesting to 
note that such scalar and pseudoscalar structures reappear in the dimension-8 effective operators, 

3 



A ^(8) = ^ E C & P > €){l>'?3 4) z p R ){mWTi L°) +H.c. 

a,/3,a',/3' 
j v 2q3 (8) 



(II.5) 

where Tj £ (1, 75), (L q ) is the left-handed lepton (quark) doublet, and <I> is the Higgs doublet 
with hypercharge | and vacuum expectation value (VEV) ($} = v/y/2(0, 1) T . We also denote, in 
Eq. |L|), = $ for g fi = d R and = ir2$* for qn = u R . Note that when the Higgs doublet 



takes its VEV, the dimension-8 operators in Eq. ( II. 5 ) reduce to the generic dimension-6 form in 
Eq. ( II. 4 ) but with coefficients further suppressed by an extra factor of v 2 /(2A 2 ). Thus, they can 
be neglected in comparison with the leading dimension-6 operators Eq. flII.4D . This means that 
for the linear realization of the SM gauge symmetry, only vector and axial-vector operators with 
r, G (7,7, 70-75) are relevant at dimension-6 level. For the phenomenological analysis in Sec. Ill, we 



will focus on analyzing the bounds for the most general form in Eq. (II.3) since Eq. ( |IL3| ) contains 



the restricted form in Eq. (II.4) as a special case. 

Finally, we note that in principle, we could also include purely leptonic dimension-6 operators 
of the form (p,Tr) (£ T£@), which contains an additional lepton bilinear (P T instead of quark 
bilinear (j[ a TqP). They may involve similar flavor dynamics as the effective operators ( |II.3j ), but 
unlike ( |II.3| ) they are relevant to only a few low energy constraints. Most nontrivial bounds come 
from certain three-body rare r decays and all appear similar. We will summarize these separately 
in Appendix B. 
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B. Theoretical Consideration for the Size of Coefficients 



The precise value of the dimensionless coefficient C^g in Eq. ( |II.3[ ) should be derived from the 
corresponding underlying theory in principle. In the current effective theory analysis, we will invoke 
a power counting estimate. As shown in the Appendix A, the operator Eq. QH.3p is formulated under 
the nonlinear realization of the SM gauge group Qsm j which provides a natural effective description 
of the strongly coupled electro- weak symmetry breaking (EWSB) sector and/or compositeness. In 
this scenario, the natural size of C^g for an effective dimension-6 four-Fermi operator such as (|II.3| ) 
can be typically estimated as ||, 

Cjg = 47r0(l), (default), (II.6) 

which corresponds to an underlying theory with a strong gauge coupling a s = <7|/47r = 0(1). 
Naive dimensional analysis (NDA) [ pl[ |22|1 provides another way to estimate operators in the 
nonlinear realization. For the dimension-6 operators in Eq. (fill), the NDA gives 



C^[NDA] < (4vr) 2 0(l). (II.7) 

which corresponds to an underlying theory with a strong gauge coupling[Q| as = <7§/47r < 
0(4ir) . So, in general, for the nonlinearly realized effective theory, we expect C J a3 > 47r, and 
to be conservative we will choose the estimate Eq. ( p.6| ) as the "default" value of our analysis. 
With Eq. ( [II, 6 ), all the phenomenological bounds derived in the next section can be translated into 



bounds on the new physics scale A. The bounds with a different counting of C J a3 , such as Eq. (|II.7| ) 
above and Eq. ( |II.8| ) below, can be directly obtained from our default results by simple rescaling. 
In the weakly coupled theories, we have 

C° aB = 0(1) , (weakly coupled scenario). (H-8) 

Since the new physics scale of a weakly coupled scenario is likely to be first determined by dis- 
covering the light new particles (such as light Higgs boson(s) and a few lower-lying states of su- 
perpartners), we may mainly motivate the current analysis by the strongly coupled theories where 
the estimate in Eq. ( |IL6j ) or ( |II.7D can sensibly apply. However, as will be clear from the Sec. Ill 
and IV below, even for the weakly coupled theories with Eq. ( P-§| ), significant bounds on the new 
physics scale A can still be derived for the operators with the quark-bilinear q a V q$ containing no 
t quarks, no c quarks and at most one b quark. 
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The linear realization is more appropriate when there is a Higgs boson with a mass well below 
the scale 4ttv ~ 3.1 TeV, such as in the typical models with supersymmetry Q or composite 
Higgs models with Top-color P||55[|. Finally, we note that the above estimates for the coefficient 
are flavor-blind, i.e., we do not worry about the possible suppression from flavor- violation effects. 
Below we will analyze how the large leptonic r — fx flavor mixing is naturally realized in typical 
scenarios without further suppression. 




FIG. 1: The charged current interactions for lepton sector involving MNS mixings (left plot), and for quark 
sector involving CKM mixings (right plot). 



C. Neutrino Oscillations, Large Lepton Mixings and r — fi Flavor Violation 
Operators 

Neutrino oscillation experiments can measure the lepton charged current interactions involving 
(large) MNS mixings (cf. left plot in Fig. 1) analogous to the quark charged current interactions 
involving (small) CKM mixings (cf. right plot in Fig. 1). Thus we can generally write, 

C cc = —% E ^ P V<xiAw- + H.c. (11.9) 

V 1 i=l,.-,N 

where V is the 3x N unitary MNS mixing matrix containing a product of two left-handed rotation 
matrices, 

3 

V ai =J2 U L/3aUZp i , (H.10) 
13=1 

where (3 x 3) is from the lepton mass diagonalization U L MfUft = M^ iag and (N x 
jV) |56|| is from diagonalizing the neutrino mass matrix (of Majorana or Dirac type) f57|| ■ Hence, the 
large or maximal mixings in the MNS matrix V can originate from (i) either UY (neutrino mass 
diagonalization), (ii) or (left-handed lepton mass diagonalization), (iii) or both sources. In the 
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cases (ii) and (in), we see that large mixings in play important roles for neutrino oscillation 
phenomena. Furthermore, the right-handed lepton rotation does not enter the MNS matrix and 
is thus free from the constraint of oscillation experiments. This means that even in the case (i) large 
lepton mixings can originate from C/jj though is constrained to have only small mixings. Clearly, 
the neutrino oscillation data alone could not identify the origin of the MNS mixings (involving 
only the product of two left-handed rotations), i.e., whether the large or maximal mixings in MNS 



matrix really originate from the neutrino mass matrix or charged lepton mass matrix or both [58]. 
Therefore, to fully understand the flavor dynamics in the lepton and neutrino sector, it is important 
to directly test lepton rotation matrices £/£ and in other lepton flavor-violation processes. The 
large or maximal MNS mixings observed in the neutrino oscillation experiments strongly motivate 
searches for large lepton flavor-violating interactions originating from and t/jj . 

The effective dimension-six t — fi operator Eq. ( pi. 3 ) can arise from exchange of certain heavy 
particles such as a heavier neutral gauge boson (Z'), or a heavier Higgs scalar (S°). The underlying 
dynamics for such interactions fall into two distinct classes: (a), flavor universal, or, (b). flavor 
non-universal. For the Class- (a), we can write a generic lepton-bilinear term tpTjip, where ifi = 
(e, (j,, t) t . We thus deduce 



Yl ^x T iA' ( r i = !> 75); 

x{¥<) 



Yl ^x T i ^x ' = 7<x , 7^75 ) ; 

X 



(11.11) 



where x, ( = L, R. The lepton mass diagonalization U^M^U^ = M^ iag enables us to transform 
leptons from flavor-eigenbasis (VO into mass-eigenbasis (ip 1 ) via 1/3 1 = Uj^ip'j^ and tpR = U^ip' R . Thus, 
Eq. ( pi.ll| ) can be rewritten as, in the mass-eigenbasis, 



ipFj if) 



£ ^ Ij (U?Ut) Vc => 0(1) t'Tj y! + 0(1) ti'Tjr', (Tj = 1, 75 ) ; 

x(^C) 

S ^'x T i ^'x ' ( r i = 7t , 7a 75 ) • 



Here, for Tj = 1, 75, the unsuppressed r' — // mixing bilinears can be induced by the large 0(1) 
"23" or "32" entries in the lepton rotation matrices U^U^ and U^U^. We see that such large 
flavor-mixings between r' and y! occur only for V = 1, 75, implying that for Class- (a) the r — \i 
operators arise from exchanging a heavy scalar (pseudoscalar) in the underlying theory (cf. Fig. 2b 
and 2c). For scalar (pseudoscalar) type couplings to be flavor universal, there should be certain 
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flavor symmetry associated with these couplings in the underlying theory. The quark bilinear 
QTjQ can join the same type of flavor universal interactions, where Q =U, T> with Li = (u, c, t) T 
and T> = (d, s, b) T . But, without fine-tuning the flavor-violations in the quark-bilinear would be 
relatively small based on the experimental knowledge about the CKM matrix, except that some 
right-handed mixings may be quite sizable. 

We then proceed to consider Class- (b) with flavor non- universal dynamics which is strongly 
motivated by the observed large mass hierarchies for leptons and quarks. For instance, the in- 
duced lepton-bilinear may contain only the third family tau-leptons as happened in the dynamical 



symmetry breaking models ||, [32], |33| and lepton non-universality models |34|, |35|, |3q| . Gener- 
ically, we can write down a tau-lepton bilinear in the flavor eigenbasis, t^Tjt^ where X = C 
corresponds to gauge boson Z' exchange and \ 7^ C corresponds to scalar S° exchange. The Z' 



exchange from the strongly interacting theories such as top-color models [32| is particularly inter 



esting for studying the effective dimension-6 operators ( pi. 3 ) as the Z' induces strong couplings for 



Eq. ( p.3| ) which naturally fit the counting in Eq. flll.6| ). Because of the flavor non-unversality of 
such underlying dynamics, the generic tau-lepton bilinear takes the following form after the lepton 
mass-diagonalization, 

T^TjT C = I^T :j (U^) 3a (U E c )^£f 0(i)7r j ii' + 0(l)'i7T j T J , (11.13) 

where x C f° r = 1> 7s an d X = C f° r Tj = r y a , 7 CT 7 5 . We see that, due to the allowed large 
entry (11^)^2 or (U^)^, the un-suppressed flavor-violating t — fi bilinear can be generated for Tj 



being either (axial-)vector or (pseudo-)scalar (cf. Fig. 2a-c), unlike the situation in Eq. (11.12). In 
most cases, such flavor non-universal dynamics generically invokes the third family quark-bilinears 
t x Tj tf and b x Tj at the same time, and some sizable right-handed mixings such as t R —c R mixing 



(or b R — s R mixing) can naturally arise [37|. To be model- independent, we will include all possible 
flavor combinations in the quark-bilinear q a Tj [cf. Eq. ( |II.3| )1 for our phenomenological analysis 
in the next section. 
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(a) (b) (c) 

FIG. 2: Illustration of typical generations of the effective four-Fermi interactions [cf. Eq. flll.3| )1 in plot (c) 
from the exchange of (a) a heavy new gauge boson Z', or, (b) a heavy new scalar or pseudoscalar S° in 
the underlying theories. Here the r' — // flavor-violation couplings solely come from the charged lepton 
mass-diagonalization matrices XJ l L and/or t/jj,. 

In summary, a full understanding of leptonic flavor dynamics (for both masses and mixings) 
requires experimental exploration not only via neutrino oscillations but also via other lepton-flavor- 
violating processes. The above classification shows that, given the large lepton mixings (especially 
for the tau and muon leptons) as motivated by the atmospheric neutrino oscillation data, the 
effective dimension-six flavor-violating t — fi operators ( |II.3| ) can be naturally realized without ad- 
ditional suppression (in contrast to the quark sector). Furthermore, particular chirality structures 
(characterized by Fj) can be singled out, depending on whether the underlying leptonic dynamics 
is universal or non-universal in the flavor space. The systematic phenomenology constraints ana- 
lyzed in the following Sec. Ill will demonstrate how the various low energy precision data on lepton 
flavor- violations can probe the new physics scale associated with the underlying flavor dynamics, 
complementary to the neutrino oscillation experiments. 

D. Radiative Corrections versus Leading Logarithmic Term 

In our analysis, we will consider two classes of bounds, from contributions at either the tree 
level or the loop level. Many of these bounds arise from the tree-level operators ( |II.3| ) directly and 
will be derived from various low energy decay channels. In some cases, the significant bounds can 
only be obtained by relating the operators involving one set of heavy quarks to those involving 
lighter quarks, through exchange of a W or charged Goldstone boson at the loop level (cf., Fig. 3 in 
Section 3). How does one handle radiative loop effects in an effective theory? For such calculations, 
typically, some loop integrals are divergent and must be cut off at a scale A. In the case of W 
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or Goldstone boson exchange, the divergence is logarithmic. We perform the loop integral by 
retaining the leading logarithmic terms in which the ultraviolet (UV) cutoff A would reliably 
represent the scale of new physics [38|. In doing the analysis, we assume only one (j2FjT)(q a Tjq^) 



operator to be non-zero[59|. Other diagrams involving closed fermion loops may have quadratic 
divergences. Unlike the logarithmic terms, such C(A 2 ) power corrections are not guaranteed to 
always represent the real contributions of the lowest heavy physical state of a mass M p h ys , so that 



to be conservative it is usually suggested 1 35 ] that one only uses the logarithmic terms (~ In A) 
computed in the effective theory, for representing the new physics contribution (~ lnM p h ys ) from 
an underlying full theory. We will take this approach for the loop analysis in Section 3, though we 
keep in mind that retaining only leading logarithms may possibly underestimate the new physics 
loop-contributions if the C(Mp hys ) terms are not vanishing in a given underlying theory. This 
exception occurs only when the heavy mass effect in the underlying theory does not obey the 
usual decoupling theorem [41||K]]. Thus, extracting the possible nonzero C(Mp hys ) terms is a 
highly model-dependent issue and is hard to generally handle in an effective theory formalism. 
The traditional "leading logarithm" approach provides a conservative estimate for the effective 
theory analysis and is justified for those underlying theories in which the effects of the heavy states 
(integrated out from the low energy spectrum) exhibit the decoupling behavior. 



III. PHENOMENOLOGICAL CONSTRAINTS 



We consider the general r — \x operators in Eq. ( |II.3| ) and take Tj to be the same for both the 
lepton and quark pieces. There are four types of operators to be considered, Tj = (S, P,V, A), 
and for each type there are twelve combinations of q a qp, (uu,dd,ss,cc,bb,tt,ds,db,sb,uc,ut,ct). 
This gives a total of 48 operators for our analysis. We first consider operators involving two 
light quarks (uu,dd,ss,ds), then the non-diagonal operators involving one or more heavy quarks 
(db, sb, uc, ut, ct), and finally consider the diagonal operators involving two heavy quarks (cc, bb, tt). 
In our analysis, we will consider one operator to be nonzero at a time and derive the corresponding 
bound on the new physics scale A. This should provide a sensible estimate of the scale A under the 
naturalness assumption mentioned in Section II, which states that there is no accidental cancellation 
among the contributions of different operators. 
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A. Operators with Two Light Quarks 



For operators with two light quarks, the neutrinoless decay of the r into a fi and one or more 
light mesons will provide the best bound. First, we establish our conventions: The PCAC condition 
for the pseudoscalar octet gives 



jj? )b (0)\nP))=iS ab ^P», (HI.14) 



where is the meson decay constant and the Cartesian components of the axial vector current 
are 

i M5)6 = 97 M 7 5 y9- (IH.15) 

Here <p = (f) a X a /V2 with a = 1 • • • 8 is the 3x3 matrix of pseudoscalar meson fields and A™ are 
the Gell-Mann matrices normalized according to Tr(A a A 6 ) = 25 ab . For current quark masses, we 
choose 

m u = m d = 5 MeV , and m s = 120 MeV , (III. 16) 

and the results are not particularly sensitive to these choices. For simplicity, the muon mass 
(and the pion mass, when applicable) will be neglected relative to the tau mass when calculating 
kinematics. 

Knowing the vacuum transition matrix element, we can readily evaluate the particle decay 
width. For instance, for a two-body decay r — > jiM where M is a generic light meson, we have the 
spin-summed and averaged partial width 

r(r - ixM) = J-J-^|M| 2 , (111.17) 

where the transition amplitude iM =< fiM\iH\r > can be evaluated by vacuum insertion. 

• Axial Vector Operators 

Bounds on the uu and dd axial operators can be obtained by looking at r —> [iir . Using 



(0\u^j 5 u\tt°(p)) = (o 



'(/>)> =»^|p„ (HI.18) 



with F n = 131 MeV, and noting that for u — > d the right-hand side is the same except with an 
opposite sign, we find 

r (r ->■ ^vr°) = —lF%m%. < °- 908 x 10_17 GeV , (111.19) 
12 



where the inequality comes from the 90%C.L. experimental bound on this decay mode listed in 
Ref. [43]. This then implies that, for both the u7 M 7 5 u and d^^d operators, 



A > 11.3 TeV. 



(111.20) 



For the ss axial operator, a bound is obtained from r — > \xr). We have 



(111.21) 



where F v = F® 7=^7 i s defined using ( 

V2 



v(p)) = i ~wP^ Here > A ° = y 1 and 



J ~^" 7 " y 3 

at the next-to-leading order (NLO) of the chiral perturbation theory [44|, F® = 154 MeV and 



F2 = 25 MeV. Using these values we find that F v ~ i^, the value from the SC/(3) limit. This 
gives, 



r(r -> /X77) 



2vr K - «^) 2 ^ 



< 2.179 x lO -17 GeV, 



A 4 



(111.22) 



implying 



A > 9.5 TeV. 



(HI.23) 



Att — 

Note that, for the isospin-invariant effective operator (ft At) (uAu + dAd) with A = 7^7 5 , the 
same bound of A > 9.5 TeV can be derived from the above process. 
For the sd operator, the bound comes from r — ► fiK . Thus, we have 



(0|s7^7 5 rf|^°> = iF K pV, 



(111.24) 



where experimentally Fk = 160 MeV. This leads to 



r(r -► llK q 



7T 



[mi — mid i* 1 ^ 1K 

K < 2.27 x 1(T 15 GeV. 



A 4 



(111.25) 



and thus, 



A > 3.6 TeV. 



(HI.26) 



• Pseudoscalar Operators 

Here, the Dirac equation is used to reduce the axial vector matrix elements to pseudoscalar 
matrix elements, and then we use the same processes as above. 
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We find that 



(0 \u^u\ n°(p)) = - (0 \d 7 5 d\ 7r°(p)> = 4= ^ F n , (111.27) 

V2 rn u + m d 



2 



which then yields 

n- ^ „-"•■ - _. _ 

8 A 4 m2 



r(r - ^vr°) = ^ F y? T < 0.908 x 10~ 17 GeV , (111.28) 



so that 

A > 11.7 TeV. (111.29) 
For the strange quark operator, we find 



2 

rat 



<0|»A| » W > = -»V6 f « mu + m « +4m< (111.30) 



which gives (taking rj = rjs) 



4 12 2\ 2 

67r fnti^n [mz — rat) 

r(r -► W ) = ^ 5-2 q „ ^ < 2.18 x HP 17 GeV , (111.31) 

^{m u + ra d + Ara s ) 2 va% 



implying 



A > 9.9 TeV. (111.32) 



47T 

For the isospin- invariant effective operator — (fiPr) (uPu + dPd) with P = 7 s , the same bound 
of A > 9.9 TeV can be derived from the above process. 
Finally, we have 

(0 \sj 5 d\ K°) = i F K , (111.33) 

rad + ra s 

which gives 

,2 



o- _ 7T m 4 K F K 2 {m 2 T - ra 2 K ) 



r(r -► jtxif ) = T77 2 v r — 3-^- < 2.27 x 10~ 15 GeV . (111.34) 



A 4 (m^ + ra t 



Thus, we deduce 

A > 3.7 TeV. (111.35) 



• Vector Operators 



We take a simple SJ7(3) relation: 



v f 



(111.36) 
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*1 

V2 



where = V^—^= is the vector meson octet. Using vector meson dominance [45|, we determine the 



dimensionless ratio g = m /c from V 



e 1 e 



for each vector V 



which yields g p = 5.1, 



guj = 17 and g^ = 12.9. These phenomenological values indicate some SJ7(3)(3) breaking, as 
expected. Assuming ideal mixing, we get 

(0\u^u\p°) = - (0\d^d\p°) = ie»K p , 

(Q\^fs\4>) = ie^K^, 
0\sYd\K *) = ie^K K , , 



where 



m, 



9p 



3mj 



and in the SU (3) limit, one has Kk* = K p . Then, we find 

.2 \ 2 



r( r - pV) = 
which gives bounds as follows 



irKy m T 



A 4 m 



nil 



r(r -> /xp) < 1.43 x 1(T 17 GeV 
T(r -> //(/.) < 1.59 x 10~ 17 GeV 
r(r -> ftK*) < 1.7 x 10~ 17 GeV 



A > 12.4 TeV, 
A > 14.3 TeV, 
A > 12.8 TeV, 



(uu, dd) , 
(ss) , 
{ds). 



(111.37) 



(111.38) 



(111.39) 



(111.40) 



• Scalar Operators 

Scalar operators will lead to three-body decays of the r into a p and two mesons. Using the 
leading order chiral Lagrangian, we obtain the matrix elements of scalar densities at the origin 

(0\ss\K + K-) = (0\uu\K+K-) = B , 

(0|nu|7rvr) = (0\dd\Tnr) = B , 
3 

3(0|mi|77 8 7? 8 ) = - (0 1 ss | r7 8 77 8 ) = B , 

(0\ds\TT + K-) = B , (111.41) 

where m 2 = 2mBQ and fn 2 K+ = (m u + m s )Bo with m = m u + m^. We take m M = = 5 MeV, 
which gives Bq = 1.96 GeV and m s = 120 MeV. Thus, the differential decay widths are computed 
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as 



1 R 2 



dT(r - ^vr^ u ) = -dT(r - = ^^-^ ^ + m^dExdE* , 



2 " 647r3A uu , d / 

rfr(r^m) = (^) 647r 3^4 ^ E 1 + m li )dE 1 dE 2 , 

dT(T - Mif+K-) = i? % 9 3^ (#i + m^)dE 1 dE 2 , 

6171 A dd,ss 

dT(T - /xK+tt-) = ^o 2 327r 3 A 4 (^i + m^)dE 1 dE 2 . (111.42) 

The best bounds on the Cm and (id operators come from 

T(r -► ^tt+tt-) < 0.186 x 10" 16 GeV, =► A uu4d > 2.6 TeV, (111.43) 

the bound on the ss operator from 

r(r -> fiK + K~) < 0.341 x 10~ 16 GeV, A uu , ss > 1.5 TeV, (111.44) 

and the bound on the sci or ds operator from 

r(r -► fiK + TT-) < 0.168 x 10~ 16 GeV, =► A ds > 2.3 TeV. (111.45) 

In summary of this subsection, we see that the bounds on operators involving two light quarks 
range from 1.5 to 14.5 TeV. Improvement in the experimental limits on the branching ratios, of 
course, will increase these bounds via the fourth root of the branching ratio. Much better im- 
provement can be obtained from processes involving decays of heavy quarks, which we consider 
below. 

B. Non-diagonal operators involving one or more heavy quark 

Now, we analyze the operators involving (uc, ut, ct, db, sb) quarks. The bounds involving an up 
quark and a charm quark are problematic since the D° can not (barely) decay into [it because of 
kinematics. The bounds on these operators will be discussed at the end of this subsection. We 
first turn to the B meson decays. 
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• B meson decays 

Using the techniques described in the above subsection, one can bound the pseudoscalar and 
axial vector operators for bd quarks by looking at B° — ► fir. For the axial vector operator, we find 
that, using the experimental limit of 8.3 x 10 -4 on the branching ratio 

r i B^,r) = I J^( 1 -^)\ - A > 8.2 TeV, (111.46) 
where we take fs = 200 MeV. For the pseudoscalar operator we find, 

nB^»r) = 8 -^(l-™P) 2 , =* A > 9.3 TeV. (111.47) 

In this latter case, we have used the result from Sher and Yuan |46|], | (0\dj 5 b\B) \ 2 = 4ms/?^/7r 3 / 2 , 
where (3b ~ 300 MeV is a variational parameter. 

In a moment, we will consider the scalar and vector operators, but let us first look at the 
pseudoscalar and axial vector cases for 6s quarks. Here, precisely the same analysis as for B — > \xr 
can be done for B s — ► /ir, with the same result for the width, in the approximation where the 
masses of the constituent s and d quarks are equal. Alas, there are no published experimental 
bounds for B s —* fir. Note that the lifetime of the B s is given by 1.46 ± 0.06 picoseconds, 
compared with the B lifetime of 1.54 ± 0.02 picoseconds. There are consistent, as expected. But, 
if the rate for B s — * fir were too large, then the lifetime would be substantially shorter. A 10 
percent branching ratio would shorten the lifetime by about 0.14 picoseconds, which would lead to 
a significant discrepancy. Without a detailed analysis, one can just conclude that there is a bound 
of five to ten percent on the branching ratio for B s — ► fj,r; we will give bounds assuming it is ten 
percent. 

With a 10% bound, the above results scale as the fourth root of the branching ratio, giving a 
bound on the axial vector operator of A > 2.5 TeV and on the pseudoscalar operator of A > 2.8 TeV. 

Note that here is a place where an experimental bound on B s — ► fir would be very useful. This 
decay is particularly important because all of the quarks involved are second and third generation, 
and new physics effects might be substantial (especially if related to symmetry breaking); this 
decay also conserves "generation" number, and is thus particularly interesting. 

We now turn to the scalar and vector operators. Here the matrix elements {K\s^b\B) and 
{Ti\d^^b\B) are needed, along with their scalar counterparts. The vector matrix elements have 



been calculated in a quark model by Isgur, Scora, Grinstein and Wise [47]. They note that for a 
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light pseudoscalar meson X, 

(X(p x )\qrb\B(p B ))^f + (q 2 )(p B +p x r + f.(q 2 )(p B -p x r, (111.48) 

and present expressions (in their Appendix B) for /+ and /_. Here q 2 = (p B — Px) 2 an d B and X 
are on-shell with X ~ qd. The masses in these expressions are constituent quark masses, which we 
take to be m q = 300 MeV, and we also take their variational parameters, 0, to be 300 MeV. For 
instance, these values give 

~m x - E x ~ 



f+(q ) ~ -/-(« ) ~ -z-J — ex P 

8 y m q 



(111.49) 



2K 2 m x 

where k ~ 0.7 is a relativistic compensation and where m x ~ 2m g . As a result of these approx- 
imations, the matrix elements should be taken cum grano salis, with an error that could be a 
factor of 2 — 4 (which translates into a factor of 1.2 — 1.4 uncertainty in the final results for A). 
The uncertainty might be somewhat larger for the matrix elements involving the pion, since the 
relativistic compensation factors are suspect. 

The result for the vector couplings is given (illustrating the B — > Kfir case) by 

dT \M 2 \ 



dE^dE-T 32TT 3 m B 



where 



\M \ = —. — - {AxA 2 - A 3 ) exp 
A 4 m d 



m B - E T - Ea- m K 



(111.50) 



(111.51) 



m d 

with A x = m 2 B -m 2 K -m 2 T - 2m B E^ A 2 = m 2 B -m 2 K + m 2 - 2m B E T , and A 3 = m 2 K [-m 2 B -m 2 T + 
m\ + 2(E IM + E T )\ . 

What are the experimental bounds for B — > Kjj,T and B — > 7T//T? None are listed. If the 
r decays semi-hadronically (which occurs 65 percent of the time) then B — » K/it will look like 



B —>■ X c \iv [50]. Then measurements of B — > X c \iv would give a higher rate than for B — * X c ev. 
These have been measured separately, with accuracies better than 0.5%, and thus an excess of 
^i-like events have not been seen with a sensitivity of 1.5% at 90% confidence level. If one assumes 
that the probability of classifying B — » K/it decays as B — > X c \iv is smaller by a factor of two, 
then, folding in the 65% branching fraction into hadrons, one would get a limit of 1.5% divided by 
0.5 for acceptance and 0.65 for the branching fraction, which is about 5 percent. A very similar 
argument would apply to B — > irfir. Obviously, a more detailed analysis could yield a substantially 
better bound. However our result only scales as the fourth root of the branching ratio bound and, 
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with the relatively large uncertainty in the matrix elements, one probably can't do much better. 
With a 5 percent branching ratio, we find that the bound on the (bs) vector operator is 2.6 TeV, 
and the bound on the (bd) vector operator is 2.2 TeV. 



For the scalar operator, one must differentiate the vector matrix elements in Eq. ( III. 48 ), taking 
care to properly include the exp [— i{ps — Pk) ■ x] factors. We find that, for instance, 

(K(p K )\sb\B(p B )) « — [f + (q 2 )(m 2 B - m\) + f-(q 2 )(m 2 B +m 2 K - 2m B E K )] , (111.52) 
nib 



where f±(q 2 ) are given in Eq. (III. 49). Using this matrix element, we find the bound on the (bs) 
scalar operator is 2.6 TeV, and similarly that for the (bd) scalar operator is 2.2 TeV. One should 
keep in mind the relatively large uncertainties in these bounds due to the hadronic uncertainties 
discussed above. 

• Top Quark decays 



Due to the very short life-time of a top quark, bound state top mesons do not exist. Bounds 
on operators involving a top quark can be readily obtained by looking for decays t — ► cprr and 
t — > ufj,r. Neglecting the final state masses, one finds that the width is 



96ttA 4 ' 
24vrA 4 ' 



for scalar and pseudo scalar couplings, 
for vector and axial vector couplings. 



(HI.53) 



CDF [48] measures the ratio, 



Br(t -> Wq) y ' 

by counting 6-tagged top events and all top to Wq events. The result is TZ = 1.23^'^, which 
translates into TZ > 0.72 at 90% C.L. For the t — * u/xt channel, one considers this to be similar to 
Wq, in that the signature is an isolated muon with some jet activity (clearly, a detailed analysis by 
CDF/DO could distinguish between Wq and ufir), and thus one approximately has Br(i — * ufir) < 
1 — 0.72 = 0.28 at one standard deviation. This leads to a constraint 



A > < 



190 GeV, for scalar and pseudoscalar couplings, 

(111.55) 

270 GeV, for vector and axial vector couplings. 



It should be kept in mind that the above bounds on A are so close to the top quark mass that 
the use of the effective field theory is not reliable. As discussed in Sec. IIB, for models obeying the 



19 



(a) 



(b) 



FIG. 3: Diagrams which relate couplings involving w-type quark bilincars to d-type quark bilincars. The ir 
is the charged Goldstonc boson. 



Naive Dimensional Analysis (NDA), the corresponding bounds become stronger than the above by 



be more reasonable. One could improve on the above bounds significantly from non-observation 
of t — > jet + /xt decay, but this has yet to be done. (We also recall that our final limits on A only 
vary as the fourth root of the branching ratio.) 

• Loop Contributions 

Operators involving heavy quarks are harder to constrain with meson decays. However, one- 
loop contributions via exchange may mediate the transition from a heavy quark to a lighter 
one. This leads to processes with external light quarks and thus results in possibly significant 
constraints from light mesons. For the vector and axial vector couplings, strong bounds can be 
obtained by considering the loop contributions with W ± exchange. Such contributions will also 
give good bounds on the (cu), (cc), (tt) operators, as will be discussed below. 

Consider the loops in Fig. 3(a) and (b), where 7r~ is the charged Goldstone boson. These 
loops will relate, for instance, a heavy-quark operator of the form Cj^/ 'A 2 (jlTkT)(Uj 1 Ti t Uj 2 ) to a 
light-quark operator of the form Cf ii2 /A 2 (j2 r)((ij 1 Tidi 2 ) . Since we have very strong bounds on 
operators with light d\ quarks, this gives a method of deriving bounds on heavy Uj quarks, and in 
many cases will provide the only bounds. 

We consider the diagrams of Fig. 3(a) and (b), and derive the quark-bilinear, 



about a factor of %/47r ~ 3.5, 



and thus in this case the application of effective theory formalism will 



5£% {qiT k qj) , 



(111.56) 
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where 5C 1 ^ is the loop-induced form factor computed from the diagrams. It is easy to show that 
there will be no contribution to the scalar and pseudoscalar operators from the loop, and thus only 
vector and axial-vector operators are relevant. 

For Tf, = 7^, labelling the two external quarks with indices i^j and the two internal quarks 
with indices (ji, j'2), we find that the vector and axial- vector couplings are generated with corre- 
sponding induced form factors, 

where Vij is the relevant CKM matrix element. We note that if a mass-independent renormalization 
procedure is used (cf. footnote in Ref. [59]), we would extract the leading logarithmic term as 
ln(A///o) with set as the energy scale of the relevant low energy process (e.g., /Uo = m& for b- 
decays). But, for deriving the final bound on A, this does not make much difference as it only goes 
like -y/ln(A/m,fc) [cf. Eq. ( [II.59Q below]. Such difference would not be a main concern for leading 
logarithmic estimates since in the leading logarithm approximation, all unknown non-logarithmic 
terms are dropped by assuming the absence of accidental cancellations. 
For Tk = 7^7 5 , we find the analogue of the above equation, 

«** = _ scr = _t^M L + _ a . (m . 58) 

v A 32vr 2 I 2M^ ) M w y ' 

We will ignore small masses of light quarks and the overall signs are also irrelevant. So, the above 
leading logarithmic contribution is universal except when the internal loop fields are both top 
quarks. It is useful that a single vector (or axial-vector) coupling induces both vector and axial- 
vector vertices. This allows us to use either the vector- or axial- vector type of light-quark bounds 
to constrain both the vector- and axial-vector type of heavy-quark operators. (Note also that in 
the limit where the muon mass is neglected, our VV and AA bounds from r decays give identical 
bounds on VA and AV.) 

In particular, the scales A associated with heavy and light operators, respectively, can then 
be related. Letting (X, Y) be the two internal quarks and (x, y) be the two external quarks, we 
derive, in a transparent notation, 



where the +(— ) sign is for the axial-vector (vector) coupling. Since, as discussed in Section 2, 
the loop cutoff A in the leading logarithmic terms can reliably represent the physical cutoff of the 
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effective theory p8| ], we may set the above logarithmic cutoff A equal to the light-quark bound 
A^ TXy . The light-quark bound A.^ TXy varies in the range around 1.5 — 15 TeV, and we may typically 
set it as 10 TeV. 

Now, we examine the vector- and axial-vector couplings for the t — c and t — u quark bilinears. 
With the internal quarks being the t and c (or u), the best bound comes from setting the external 
quarks to be b (which attaches to the internal t line) and d (which attaches to the internal com 
line). We then use the strong axial- vector bound 8.2 TeV for b — d quark bilinear, to obtain the 
bounds 310 GeV for the t — c operator and 650 GeV for the t — u operator which hold for both 
axial- vector and vector type of couplings. This bound is much stronger than that from the top 
decays. 

• Charm Quark off-diagonal operators 

The relevant charm quark operator is the c — u operator. One can obtain bounds from the loop 
corrections discussed above for the vector and axial-vector operators. There are several possible 
choices for the external lines (and corresponding ss- vector, dd- vector or sci-axial- vector bounds); 
the best comes from the vector ss-bound, i.e., from r — > n<j) . We find that the bound on both the 
vector and axial-vector operators is 550 GeV. 

The loop corrections will not give a significant bound on the scalar or pseudoscalar operators 
due to the chirality structure of the couplings. We do get contributions from the finite parts of the 
loop integrals, but the bounds are well below the W mass and thus not useful. For the pseudoscalar 
operator, one would get a strong bound from D° — > fir, if it were kinematically accessible, but 
it falls 20 MeV short. We have also considered virtual r decays in D — > ^r* — ► ^v+hadron, 
but this width is proportional to T T so that the realistically obtainable experimental bound on 
D — > z/r+hadron will not give anything more than a few GeV bound for A. The scalar operator 
would require an additional pion in the final state, which makes it even more inaccessible. We 
know of no bounds on the scalar and pseudoscalar operators. 

C. Diagonal Operators 

• c — c Operators 

The c — c — — t vector and axial vector operators can be bounded by the loop contributions 
(cf. Fig. 3) as discussed in the last subsection, in which the internal quarks are both charm quarks. 



22 



The bound derived from r — > fj,(f> decays is A < 1.1 TeV, for both the vector and axial vector 
operators. We note that there is no experimental bound available yet for the decay J/^f — > jj,T . 

One could bound the scalar and pseudoscalar operators by looking at [it final state of Xc an d 
r) c decays, respectively. However, no experimental bounds on these decays are available yet. 

• b — b Operators 

The obvious systems to look at are the bb bound states. No experimental bound on T —* fir 
has been published. The ratio of the decay T(ls) — * [it through the vector operator to the decay 
T(ls) — ► [i + [i~ is independent of the matrix element (if the mass difference between [i and t is 
neglected) and is given by 1447r 2 M-£/(e 4 A 4 ). The T(ls) — > [i + [i~ branching ratio is 2.5% [49|. The 



upper bound on T(ls) — ► [it can be estimated by using Ref. |^9| which measured T — * tt, and 
by comparing with the measurement of T(ls) — * [i + [i~ . If one assumes universality, these will be 
equal. One can see that the excess of tt events must be less than about 0.40% at 95% C.L. One 
then asks what fraction of \ir events would pass the cuts of the tt analysis. The tt analysis selects 
events with one r decaying to evv and the other decaying to one-prong non-electron final states; 
this would be satisfied by [xt events depending on the cut on the momentum of the non-e track. A 
conservative estimate |5(J gives an upper bound of 4% on the branching ratio, so that we arrive at 
A > 180 GeV. We know of no bound on the scalar, pseudoscalar, or axial vector couplings, since 
there is very little data on these bound states. 

The loop contributions are negligible in this case, primarily because the CKM matrix elements 
for b to c, u transitions are small. 

• t — t Operators 



The only possible way to bound the t — t — [i — T operators is through the loop discussed above, 
turning the top quarks into b, s or d. The best bound comes from the case in which the external 
quarks are b and d, leading to B° — > fir. Due to the small CKM matrix elements such as Vtd, the 
bounds are not very strong. For the axial vector operator, we find that A > 115 GeV, and for the 
vector operator, A > 75 GeV. These bounds are below the mass of top quark and thus the effective 
theory description is no longer valid. As discussed below Eq. (III. 55) for the constraints from top 
quark decays, the NDA analysis does increase the bounds by about a factor of \/4vr ~ 3.5, but even 
in this case the effective theory formalism may not be so reliable. Nevertheless, these weak bounds 
(if they might be meaningful at all) are the best ones which we could obtain at the present. For 
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the scalar and pseudoscalar operators, we know of no reasonable bounds at all. 
• Radiative r decays 

One might expect that strong bounds on diagonal operators could be obtained from r — > ^7, 
where the two quarks come together to form a loop. If the photon is attached to the quark loop, 
the result for an on-shell photon vanishes. But, if the photon is attached to the tau or muon line, 
the quark loop is quadratically divergent and independent of external momentum. As discussed 
in Sec. IID, however, quadratically divergent corrections are not guaranteed to represent the real 
contributions of a heavy physical state, and may be absorbed via renormalization, and thus will 
not be considered further. 



IV. SUMMARY 

The strong experimental evidence for large Vn — v T oscillation motivates us to explore the 
allowed mixings between the second and third generations in the charged lepton sector. We have 
systematically analyzed bounds on the generic dimension-6 flavor-violation r — fi operators of the 
form 

(/zivxriY), 

where the Dirac matrices T £ (S, P, V, A) are the same for both the t — \i and the q a — q@ 
bilinears. Such effective operators are interesting as they can naturally arise from various new 
physics scenarios and reflect the underlying flavor-mixing dynamics which may directly link to the 
large Va — v T neutrino oscillation. Since the neutrino oscillations measure the MNS mixing matrix 
which is only a product of two rotation matrices from the lepton and neutrino mass-diagonalizations 
[cf. Eq. ( II. ICQ ], it is thus important to fully understand the flavor dynamics and the origin of the 



neutrino oscillation phenomena by directly testing possible large lepton mixings from additional 
flavor-violation processes. Given such generic dimension-six r — fi effective operators, we have 
considered all possible flavor combinations in the quark-bilinear q^T q@ and analyzed existing 
experimental data for a variety of processes to establish the best available bounds on the invoked 
new physics scale A. Our results are summarized in Table |[ 

In Section IIIA, we studied the operators in Eq. ( II.3[) involving two light quarks. We found 



quite strong bounds in this case since there are good experimental limits on decays of r to \x and 
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Bound 


1 7 5 7a 7a7 5 


uu 


2.6 TeV 12 TeV 12 TeV 11 TeV 

(r — > p7r + 7r~) (r — > fin ) (r — > pp) (r — > p7T°) 


dd 


2 6 TeV 12 TeV 12 TeV 1 1 TeV 

(r — > p7r + 7r~) (r — > p7r°) (r — > pp) (r — > p7T°) 


ss 


1 5 TeV 9 9 TeV 14 TeV 9 5 TeV 
(r — > fiK + K—) (r — > pr/) (r — > p</>) (r — > pry) 


sd 


2.3 TeV 3.7 TeV 13 TeV 3.6 TeV 
(r -> fiK+n-) (r -> pK°) (r -> pK*) (r -> pK°) 


Id 


2.2 TeV 9.3 TeV 2.2 TeV 8.2 TeV 
(B -> ^pr) (B -> pr) (5 -> Trpr) (5 -> pr) 


is 


2 6 TeV 2 8 TeV 2 6 TeV 2 5 TeV 
(B^Tvpr) (B s ^pr) (B -> Xpr) (B s -> pr) 


tc 


190 GeV 190 GeV 310 GeV 310 GeV 

UU V.I V V 1JU V.' V V U IV V_< V V U 1U v_J Vy V 

(t — > cpr) (t — > cpr) (B — > pr) (B — > pr) 


iu 


190 GeV 190 GeV 650 GeV 650 GeV 
(i — > -upr) (i — > -upr) (B — > pr) (B — > pr) 


cu 


* * 550 GeV 550 GeV 

(r -> p0) (r -> p0) 


cc 


* * 1.1 TeV 1.1 TeV 

(r -> p0) (r -> p0) 


bb 


★ * 180 GeV * 

(T -> pr) 


it 


* * 75 GeV 120 GeV 

(B - pr) (B - pr) 



TABLE I: Bounds at 90% C.L. on four-Fermi flavor- violation operators of the form (p,TjT)(q a Tj q 13 ), where 
Tj e (S, P, V, A). Combinations for which no bound has been found are marked with an asterisk, otherwise 
we list the process which gives the strongest bound (cf. text for details). 

one or two light mesons. For P, V and A the bounds on A are of order 10 TeV (except for sd where 
the P and A bounds are weaker since the experimental constraint is weaker). Since three-body 
decays are used for the S case, the bounds obtained are slightly smaller, of order 1.5 — 2.5 TeV. 
In Sections IIIB and IIIC, we studied operators involving at least one heavy quark. For the 
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case of a charm quark in Eq. (II. 3), one might think of using D-decays to final states involving 
fjiT but these are ruled out by kinematics. However, it turns out that we obtained good bounds 
on operators involving the cu and cc combinations for the V and A cases by considering loop 
contributions (shown in Fig. 3) to r — > fMp. Loop contributions involving c quarks are enhanced by 
the fact that V cs ~ 1. For the S and P cases we could not find a bound on A since the loop diagrams 
shown in Fig. 3 have vanishing leading logarithmic term and the remaining finite loop terms are 
numerically negligible. Also there is no bound from pseudoscalar and scalar charmonium decays 
- these states are rather broad with several MeV uncertainties in their widths at the present, so 
there is no significant experimental constraint on their branching fractions to \ir. 

For operators involving one 6-quark, we used the experimental limits on B — ► /xr, B s — > /it and 
estimated bounds on8-> iitM |6l| to obtain bounds on A in the 3 — 9 TeV range. There is some 
uncertainty in our approximate HQET|32| hadronic matrix elements which could be improved, but 
since it is a fourth root which appears in our extraction of A, the final error is not so large. For 
the bb case, we considered the contribution of loop diagrams in Fig. 3 with internal b quarks to 
processes involving external u quarks. However, these are suppressed by a CKM factor |T^h| 2 and 
thus do not lead to a useful bound. The only significant bound for bb which we could obtain is for 
the y-type operator and it comes from T decay. 

For operators involving one top quark (via tu and tc), the best bounds for the S and P cases 
come directly from the top quark decays t — ► c[at,u[j,t. Their decay widths scale as mj 5 /A 4 so 
that the bounds are actually non-trivial. For the V and A types of operators with quark bilinears 
tu, tc and also it, the strongest bound comes from the associated contribution to B — > fj,r due to 
internal t quarks in Fig. 3. We have not found a way to bound the S and P operators for the case 
with ft quark bilinear. 

It is important to note that under the linear realization of the Standard Model gauge group, 
the allowed operators are restricted to V and A types, as discussed in Sec. IIA and Appendix A. 
In this case, we obtain bounds for all but one of the allowed operators (i.e., except T = A for the 
quark-bilinear b T b) . We should mention that if the coefficients C^g of our dimension-6 operators 
follow the estimate of the NDA analysis in Eq. ( |II.7| ) for certain class of strongly coupled theories 
[instead of the "default" estimate in Eq. flH.7p ], the final bounds in Table I would be stronger by 
a factor of < y4w ~ 3.5. On the other hand, if the underlying theory is weakly coupled (such as 



SUSY-type models), the coefficients C J aj3 ~ 0(1) [cf. Eq. ( II. 8 )1 and thus the final bounds in Table I 
would be weaker by about a factor of v4tt- Impressively, even in this weakly coupled scenario, 
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Table I shows that significant bounds of A > 0.6 — 4TeV still hold for all those quark bilinears with 
no t or c quark and at most one b quark. Also, it will be interesting to further investigate how the 
present bounds in Table I can explicitly constrain the relevant models classified in Sec. IIC, which 
will induce specific forms of our fi — T flavor violation operators in the low energy theory. 

We note that most bounds listed in Table I are from rare decays of r and B. Tighter experimental 
bounds on r decays would lead to stronger bounds on almost one half of the operators considered 
here. Searches for fj,r from charmonium decays (J/^,Xc,i]c) would be an important addition to 
study this class of operators. Also, as emphasized in Section III, it would be particularly significant 
to have an experimental bound on B s — > \xt and also B decays to (itM. Since we have not found any 
way to bound three of the four operators with bb bilinears, it would be helpful to have experimental 
bounds on scalar and pseudoscalar bb decay to /ur. It is also very interesting to search for the decay 
t — > [it + jet at the top-quark factories such as the Tevatron Run-II and the CERN LHC. 

Finally, in Appendix B, we also considered an extension of our formalism in Sec. IIA to include 
purely leptonic r — fj, operators at dimension-6 with a lepton-bilinear {£ T £@) [cf. Eq. ( |B.l| )]. 
Among the available constraints, we found that the three-body rare decays r — > 3/i, ////e, \iee give 
the best bounds at the order of 10 TeV or so. For (frf) being a neutrino pair, the bounds from 
r — > jiuv decay are much weaker, around 2 — 3 TeV. No significant bound is obtained for (jf F £@) 
containing one or two heavy r leptons. 
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APPENDIX A: NONLINEAR VERSUS LINEAR REALIZATION OF THE ELEC- 
TROWEAK GAUGE SYMMETRY 



We note that in principle the exact form of Eq. (II.3) depends on how the electroweak gauge 
symmetry is realized in £sm- The general form ( |II.3D can be derived by using the nonlinear 
realization of the SM gauge symmetry. Under the nonlinear realization, the SM Higgs-Goldstone 
fields are parameterized as 

$ = _ ( V + H)U, U = exp[ivr a r a /w] , (A.l) 
a/2 

which transforms, under Qsm = SU(2)l f8> U(1)y, as 

U^U' = g L Ug Y , H^H' = H, 

g L = e W [-iO a L T a /2] , g Y = exp[-i^r 3 /2] . 



(A.2) 



We introduce the following useful notations, 



r 3 



(A.3) 



D»U = d^U + igW«—U - ig'UB^- , 

= —Tr^rfD „U) = --Tr^vg = --TrpT^] , 
9 9 9 

z o = _!^ Tr[T 3^ D U] = _!^ Xr[r 3y ] = J^L^y ] 
^9 9 9 

where c w = cos 9w and s w = sin 6\y . It can be proven that = 
[W^"t~ + VV^7 r+ + c w 1 ^ r3 ]- Then, we can write the following transformation laws, under 
Qsm, 

v„ - K = 9 L vA , - % = gyV^y , 

(A.4) 

Wf - K' = exp[=Fi0y]W£ , Zl - 2* = Z° . 
Since (W^, Z°) feel only the unbroken U(l) em gauge interaction, its covariant derivative is 

D U W± = {d u ± ieA v )W± , = d v Z° , (A.5) 

In general, the non-linear composite fields (W^, Z°) can be expanded as 

^ = + ^ V + £ Wit" - W^) + ^ - 

(A.6) 
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so that in the unitary gauge, (W^ 1 , Z®) = (W^ 1 , Z®), and 



V2W j 



(A.7) 



Now, we Cciii rewrite the SM Lcigrangian in terms of fields which feel only the unbroken £^(l)em; 
£SM = ^SM + ^SM + ^SM 

£sm = faHWH - V SM (H) - ±Tr [V*%] [v 2 + 2vH + H 2 ] (A.8) 
^sm = -\f~l^ {K tT + C ^ Z > 3 ) FjL + fr-f fa + ieQf.Al - ig's w Q fj Z°] ft 
-TjMljfil + H/v) 

where is the general fermion mass-matrix and after diagonalization, M^, = TnfjSjji . The 

electric charge of the fermion fj is defined by Q f. = I^j + Y/./2 , and Fjt = (fijL, fijL) T is the 
left-handed SU(2)l fermion doublet. For the electroweak interactions in this non-linear realization, 
Eqs. ( |A,8| ) and ( |A.4| ) show that the fermions have the U(l) em covariant derivative 

D lx f j = (d lx + ieQ fj A lx )f j . (A.9) 

Since the fermion fields only feel an unbroken electromagnetic U(l) cm gauge symmetry, we see 
that under this formalism the dimension-six t — fx operator contained in AC indeed takes the most 
general form as in Eq. ( |II.3| ) that involves one t — fi bilinear and one quark bilinear. This nonlinear 
formalism is particularly motivated when the Higgs sector of £sm is strongly coupled or the Higgs 
boson does not exist p, |12|]. In this case, the electroweak symmetry breaking scale is bounded 



from the above, i.e., A E w ^ 4nv [y, |2|. 

If the Higgs boson H is relatively light, we may choose the linear realization for the £sm [15], 



in which we consider A > Aew s ° that at the new physics scale A, the effective Lagrangian C e s 
should be invariant under SM gauge group Gsm = SU (2)l <8> U (l)y- The SM fermion fields in these 
two formalisms are related 1H 



(A.10) 

fj -» fj = exp[i9yQf.]fj , (under £ S m) , 
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where the linearly realized fermions have the usual covariant derivatives as in the SM, 

(A.ll) 

This further restricts the form of the dimension-6 operator in Eq. ( |II.3| ) because of the requirement 
of both the isospin and hypercharge conservations. Thus, we arrive at Eqs. ( II.4 ) and ( |II.5| ) for the 
linearly realized effective r — \x operators in Sec. IIA. 

APPENDIX B: BOUNDS ON THE PURELY LEPTONIC OPERATORS 

In principle, we can extend the formalism in Sec. II to include the purely leptonic operators of 
the form, 

which contains an additional lepton bilinear (£ a T £@) instead of quark bilinear (q a T q@) in Eq. ( ]II.3D . 

For both £ a and £ b being light leptons (electrons or muons), the three-body rare r decay, 
r~ — > , provides the best bounds. There are in total four allowed combinations from 

Eq. (|B.l]), (£ a , £ ) = , e + ), e~), (e~, e + ), giving the decay channels r~ — > 



(B.2) 



fi fx fi + , fi fi e + , fx /U + e , fx e e + . Their partial decay widths are computed as 

' 1, (Tj = S,P), 
^ 4, (Tj = V, A), 

where g = 1 for r~ — ► /j,~n~e + , and g> = 1/2 for t~ — > (j,~fj, + e~, fj>~e~e + . 

The experimental bounds on the decay branching ratios of these channels are given, at 90% C.L. 



^ " 768^ X < 



Br[r _ -► < 1.9 x 10~ 6 , Br[r _ -> M~//~e+] < 1.5 x 10~ 6 , 

Br[r~ -» A»~A» + e~] < 1-8 x 10~ 6 , Br[r~ -+ ^ _ e - e + ] < 1.7 x 10~ 6 , 
which are all at 10 -6 level. From these we derive the following bounds on the scale A, 

(12.8, 13.5, 11.9, 11.0) TeV, (Tj = S, P), 
(18.0, 19.1, 15.4, 15.6) TeV, (I 1 ,- = V, A), 



(B.3) 



A > < 



(B.4) 
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for the four rare decay channels, r _ — > p~p~p + , fi~p~e + , p~p + e~, p~e~e + , respectively. We see 
that these are quite similar, around 11 — 14TeV for r, = S, P and 15 — 19TeV for r, = V, A. 

For the case where £ a and are neutrinos, (£ a , £P) = {v a , v@), then the decay rate for r — > \ivv 
will be increased. Given the current data p3|| , Br[r — ► //z/z/] = (17.37 ± 0.07)%, we see that an 
increase of this branching ratio by an amount of 0.115% = 1.15 x 10~ 3 can be tolerated at 90% C.L. 
Similar to Eq. (|B.4| ), we derive the following bound from the r — ► pu a u^ channel, 

(2.2TeV, (Tj = S, P), 
A > <^ (B.5) 
[3.1TeV, (T j = V,A), 

which are weaker than ( |B.4[ ) by about a factor of (1000) 1 / 4 ~ 5 — 6 due to the different branching 
ratio. 

Finally, one or both of (£ a , £&) can be the r lepton. In this case, we may use the triangle VF-loop 



to relate the heavy r lepton operators to the neutrino operators, similar to Fig. 3 and Eq. ( 1II.59|) . 
The resulting bounds are around 150 GeV for the S and P operators, and around 200 GeV for the 
V and A operators, which are rather weak and less useful. 
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